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ABSTRACT 


New  data  from  the  Seismological  Research  Observatory 
(SRO)  is  used  in  conjunction  with  a  non-linear  least 
squares  technique  to  invert  surface  wave  group  velocity 
data  for  the  shear  velocity  structure  of  the  South-China 
subplate.  A  multi-filter  analysis  was  used  to  produce 
group  velocities  for  the  fundamental  mode  Love  and 
Rayleigh  wave.  A  number  of  earthquakes  over  a  single  path 
are  used  in  order  to  deduce  a  measure  of  the  observational 
uncertainty  of  the  present  method.  Group  velocity  standard 
deviations  range  from  .03  km/sec  to  .15  km/sec  for  Rayleigh 
waves  and  .03  to  .1  km/sec  for  Love  waves  over  the  period 
range  of  10  to  65  seconds. 

Results  of  the  inversion  indicate  that  a  3  layer, 

40  km  thick  crust  is  an  adequate  model  consistent  with  the 
data  for  the  South-China  subplate.  Crustal  shear  veloci¬ 
ties  are  found  to  be  higher  and  upper  mantle  shear  veloci¬ 
ties  lower  than  previously  reported  values  for  this 
region. 


INTRODUCTION 


The  China  mainland  consists  of  a  number  of  tectonic 
subplates,  each  of  geologically  and  geophysically  distinct 
structures.  Based  on  an  investigation  of  the  surface 
geology.  Sun  and  Teng  (1977)  divide  China  into  4  basic 
tectonic  units  (Figure  1) :  (1)  the  Ching- Tibet  subplate, 

(2)  the  North-China  subplate,  (3)  the  Northwest-China 
subplate,  and  (4)  the  South-China  subplate.  The  present 
study  focuses  its  attention  on  the  South-China  subplate. 
Specifically  it  makes  use  of  a  group  of  closely  spaced 
earthquakes  which  give  a  number  of  pure-path  surface  wave 
trains  recorded  at  a  single  SRO  station.  The  resulting 
dispersion  curves  give  a  measure  of  the  observational  un¬ 
certainty  of  the  surface  wave  data. 

Using  a  non-linear  least  squares  procedure  based  on 
the  stochastic  model  (Franklin,  1970,  Jordan  and  Franklin, 
1970,  and  Derr  et  al . ,  1970),  intermediate-  and  long- 
period  group  velocity  data  were  inverted  to  obtain  a 
detailed  shear  velocity  model  for  the  crustal  and  upper 
mantle  structure  of  the  South-China  subplate  (SOCHSP) . 
Recently  a  number  of  large  earthquakes  and  aftershocks 
have  occurred  at  the  boundary  between  the  Ching-Tibet 
subplate  and  the  South-China  subplate.  These  events  have 
been  recorded  at  an  SRO  station  at  Taipei,  Taiwan  with 
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FIGURE  I.  TECTONIC  SUBPLATES  OF  CHINA 


the  station  code  TATO.  SOCHSP  is  ideally  suited  for  a 
surface  wave  study  for  three  reasons:  (1)  the  South- 
China  subplate  can  be  considered  nearly  laterally  homo¬ 
geneous  (Tung/  1974;  Sun  and  Teng,  1977);  epicentral 
distances  are  adequate  to  allow  an  intermediate  to  long 
period  dispersive  wave  train  to  develop;  (3)  .  *ismic 
events  occur  and  are  recorded  within  the  subplate,  thus 
eliminating  the  effects  of  crossing  plate  boundaries. 

Lacking  good  geophysical  data,  most  of  the  previous 
seismological  studies  were  either  limited  to  a  simple 
estimate  on  the  crustal  thicknesses  and  shallow  crustal 
structures  or  lacked  a  discussion  of  model  resolution  and 
standard  deviation.  Knowledge  about  the  crustal  and  mantle 
shear  velocity  structure  of  SOCHSP  is  poor.  Tseng  and  Sung 
(1963)  determined  crustal  thickness  by  applying  the  tri¬ 
partite  method  (Press,  1956)  to  two  events  from  the  New 
Britain  Islands  for  phase  velocity  information.  The 
inversion  results  of  their  study  are  seen  in  Figure  2. 

The  numbers  inside  each  triangle  represent  the  average 
crustal  thickness  within  three  stations.  Using  short 
period  (4-12  seconds)  Rayleigh  and  Love  wave  group 
velocity  dispersion,  Sung  et  al.  (1965)  determined  the 
thickness  of  the  sedimentary  layers  for  several  regions 
in  China.  Tung  (1974)  made  use  of  long  period  group 
velocity  data  and  tried  to  obtain  a  shear  velocity 
structure  for  SOCHSP.  But  the  method  he  used  was 
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essentially  a  trial-ana-error  surface  wave  inversion 
technique  and  thus  he  was  unable  to  give  a  discussion 
about  resolution  or  standard  deviation  of  his  resulting 
models. 

Previous  work  has  illustrated  that  Rayleigh  waves 
produce  more  information  about  earth  models  than  Love 
waves  (Wiggin,  1972;  Derr  and  Landisman,  1972).  Wiggins 
(1972)  has  shown  that  by  simultaneously  inverting  Love 
and  Rayleigh  wave  data,  it  is  possible  to  obtain  both 
shear  wave  and  density  structures.  Though  group  velocity 
and  phase  velocity  contribute  almost  the  same  information 
(Bloch,  1969,  Bloch  et  al . ,  1969,  Wiggins,  1972),  Bloch 
(1969)  and  Fix  (1975)  come  to  the  conclusion  that  group 
velocity  is  more  sensitive  to  the  earth  model  and  there¬ 
fore  should  produce  better  resolution  in  the  inversion. 
Through  a  multiple  filtering  technique,  this  study  first 
obtains  group  velocity  measurements  from  the  fundamental 
mode  Love  and  Rayleigh  waves,  followed  by  a  Gilbert- 
Backus  (1967,  1968)  type  inversion  to  obtain  a  crustal 
and  upper  mantle  structure  of  SOCHSP  together  ^ith  an 
estimate  of  model  standard  deviation  and  resolution. 
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SOUTH-CHINA  SUBPLATE 


The  South-China  subplate  (SCS)  is  bounded  on  the 
south  and  east  by  the  Pacific  Ocean,  on  the  west  by 
longitude  103 °E,  and  on  the  north  by  latitude  32 °N 
(Figure  3) .  From  previous  gravity  and  seismic  studies 
(Tseng  and  Sun,  1963;  Tseng,  1973;  Tung,  1974)  the  crustal 
thickness  of  SCS  is  believed  to  increase  from  30  km  at 
the  eastern  boundary  to  50  km  at  the  western  boundary  with 
an  average  thickness  of  40  km. 

The  subplate  is  composed  of  two  parts,  the  Yangtze 
fault  block  in  the  northwest  and  the  South-China  folding 
block  in  the  southeast.  The  Yangtze  fault  block  is  a 
section  of  old  continental  crust  formed  during  the  Jinning 
Orogeny  (1300-950  m.y.),  and  was  subsequently  faulted  and 
strongly  metamorphosed.  Formed  during  the  Caledonian 
Orogeny  (570-400  m.y.),  the  South-China  folding  block 
was  a  section  of  transitional  (miogeosynclinal)  crust 
that  later  underwent  weak  metamorphism  and  folding. 
(Tectonic  Map  Compiling  Group,  1974,  Guizhou  Geologic 
Team  103,  1975) . 

While  surrounded  by  a  number  of  seismically  active 
zones,  SOCHSP  is  one  of  the  most  stable  regions  in  China. 
Though  a  number  of  large  events  occur  on  or  beyond  its 
western  boundary,  which  marks  the  India-Eurasia  collision. 
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south-china 


there  are  only  a  few  interplate  earthquakes  (Figure  4) . 
According  to  Li  et  aJL.  (1974)  ,  this  can  be  explained  by 
the  fact  that  this  region  has  undergone  weaker  metamorphisra 
than  other  regions  and  therefore,  the  rocks  are  softer  and 
do  not  lend  themselves  to  high  stress  accumulation. 


DATA 

The  data  used  in  this  study  are  surface  waves  from 
four  earthquakes  (Table  1)  focused  in  the  Szechwan 
Province  of  mainland  China.  The  events  were  recorded  at 
SRO  station  TATO.  A  number  of  events  over  one  path  were 
desired  in  order  to  establish  data  repeatability  and  to 
obtain  observational  uncertainties.  Because  of  instrument 
problems  and  poor  dispersion  results  (due  to  the  arrival 
of  wave  energy  at  certain  periods  that  have  taken  non¬ 
least-time  paths) ,  certain  components  of  individual 
earthquakes  were  not  used  in  obtaining  group  velocity 
information.  The  epicentral  region  and  path  are  shown 
in  Figure  3  (4-4'). 

SRO  data  in  digital  form  on  tape  is  convenient  for 
computer  processing  techniques  and  the  storage  of  a  large 
amount  of  information  on  a  single  data  tape.  The  SRO 
station  operates  at  high  sensitivity  thus  expanding  the 
detection  threshold  for  seismic  events  making  small 
magnitude  events  (M*v5)  useful  for  surface  wave  analyses. 
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TABLE  1 


LIST 

OP  SEISMIC 

EVENTS 

Date 

Location 

Focal 

Depth 

(km) 

Magnitude 

Epi central 
Distance 
(km) 

August  19,  1976 

32.893N 

33 

5.4 

1902.0 

104.189E 

August  22,  1976 

32. 992M 

33 

4.7 

1884.0 

104.181E 

September  1,  1976 

32.460N 

18 

5.1 

1886.0 

104.152E 

September  3,  1976 

28.040N 

33 

5.2 

1854.0 

100.345E 


With  this  new  system,  the  amount  of  data  available  for 
long-period  surface  waves  studies  is  greatly  increased. 

SRO  datatapes  contain  data  from  three  long-period 
channels,  sampled  at  the  rate  of  once  per  second  and  one 
vertical  component  sampled  at  the  rate  of  twenty  times 
per  second.  Long-period  data  is  recorded  continuously 
on  tape;  short-period  data  is  recorded  only  when  an  event 
has  been  detected,  and  is  sampled  at  the  rate  of  twenty 
times  a  second.  The  vertical,  north-south,  and  east-west 
long  period  channels  are  multiplexed  into  each  record. 

Each  SRO  data  record  consists  of  1000  words  (2000  bytes) . 
The  first  ten  words  are  header  information.  The  next  990 
words  are  data;  therefore,  there  are  5  minutes,  30  seconds 
of  data  per  long  period  record  (Peterson  et  al. ,  1976). 

The  amplitude  and  phase  response  for  the  long-period 
instrument  is  given  in  Table  2.  An  analytic  form  of  the 
transfer  function  used  in  our  computer  program  is  found 
by  numerically  fitting  the  amplitude  and  phase  response 
into  a  polynomial.  The  transfer  function  has  the  form; 

j  „i 

b  Cl  •  ii 

i  1 

n  i 

l  b.ZA 

i  1 

where:  a0  =  1.7443  x  10" H 
a2  =  .21163 


at  =  3.1635  x  10"3 
a  3  =  .14127 
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0667 


.05556 

.04546 

.04 

.0333 

.02778 

.02326 

.02 

.01667 

.0125 

.01 

.00667 

.005 

.0025 

.00167 


.460 

.672 

.912 

1.00 

1.01 

.9 

.724 

.550 

.400 

.200 

.108 

.031 

.0104 

.00084 

.000184 


-1.33 

-1.00 

-0.63 

-0.42 

-0.13 

0.10 

0.30 

0.47 

0.63 

0.92 

1.10 

1.34 

1.51 

1.84 

2.00 
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b0  =  3.531  x  10“  bt  =  3.3316  x  10~2 

b2  =  .12378  b3  =  1.0 

A  plot  of  this  transfer  function  is  given  in  Figure  5. 

Each  seismic  event  must  be  decoded  and  plotted  from 
the  data  tapes.  Pure  Love  waves  are  obtained  by  rotating 
the  data  coordinates.  An  example  of  a  plotted  seismic 
event  before  and  after  coordinate  rotation  is  given  in 
Figure  6  and  7. 

A  multiple  filtering  analysis  (Dziewonski  and  Hales, 
1970;  Hermann,  1973;  Tung,  1974;  Seekins  and  Teng,  1976) 
has  been  applied  to  the  surface  wave  train  to  obtain 
group  velocity  dispersion  curves.  A  typical  contoured 
result  of  this  type  of  analysis  is  shown  in  Figure  8. 
Because  of  the  epicentral  distances  involved,  only  periods 
of  up  to  65  seconds  were  used  for  Rayleigh  waves  and  50 
seconds  for  Love  waves.  The  dispersion  curves  for  the 
fundamental  mode  Love  and  Rayleigh  waves  are  plotted 
in  Figures  9  and  10.  In  these  figures,  the  solid  line 
represents  the  average  group  velocities  obtained  from 
the  multiple  filtering  analysis.  The  statistical 
properties  of  the  observations  are  listed  in  Table  3  and  4. 

THEORY 

Generalized  least  squares  methods  for  determining  the 
Earth's  structure  from  seismic  data  have  been  formulated  by 


13 


FREQUENCY 

FIGURE  5 


CFCn*OS 

FIGURE  7 

ROTATED  SEISMOGRAM  8/18/76 


TABLE  3 


RAYLEIGH  WAVE 

AVERAGED  OBSERVED  ESTIMATED 


PERIOD 

GROUP  VELOCITY 

64.00 

3.61 

51.20 

3.54 

42.667 

3.44 

36.571 

3.35 

32.00 

3.23 

28.444 

3.14 

25.60 

3.07 

21.333 

2.95 

19.692 

2.91 

13 . 236 

2.06 

17.067 

2.83 

16.00 

2.77 

15.059 

2.78 

13.474 

2.83 

12.19 

2.83 

11.13 

2.83 

10.24 

2.31 

(KM/SEC)  STANDARD  DEVIATION 

.15 

.07 

.04 

.03 

.03 

.03 

.03 

.06 

.06 

.05 

.05 

.07 

.07 

.03 

.06 

.06 

.03 
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TABLE  4 


PERIOD 

51.2 

42.667 

36.571 

32.00 

28.444 

25.60 

23.273 

31.333 

19.692 

18.286 

17.067 

16.00 

15.059 

13.474 

12.19 

11.13 

10.24 


LOVE  WAVE 


AVERAGED  OBSERVED 

GROUP  VELOCITY  (KM/SEC)  STANDARD  DEVIATION 


3.19 

3.08 

3.09 


.04 

.06 

.03 
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numerous  authors  (Jackson,  1972;  Wiggins,  1972;  and 
Crosson,  1976).  The  basic  logic  in  the  inversion  theory 
is  as  follows; 

The  observations,  Y^,  are  related  to  certain  model 

parameters,  X.,  in  some  known  way,  Y.  =  A(Xi ,X2 , . . .X  ) . 

3  i  P 

A  quasi-linear  relationship  is  established  by  taking  a 
Taylor  series  expansion  about  some  initial  value  XJ . 

»  3A,  , 

Y.  =  A. (X„)  +  _ ±  V«AX  +  higher  order  terms 

1  10  3Xj  x° 

Y±  H  Ai(x;j  +  AYi 

Ignoring  the  higher  order  terms,  we  obtain  the 
functional  relationship; 

3Aj[ 


iYi  =  ts^  ix;«j 


(i) 


where  AY^  =  observed  data  minus  calculated  observations, 
Equation  1  in  matrix  form  becomes; 


AAX  =  AY 

Where  A  is  a  matrix  whose  elements  are  the  partial 
derivatives  in  Equation  1,  AX  is  a  p  x  1  vector  whose 
elements  are  correction  to  the  initial  model  vector  xj, 
and  AY  is  the  difference  between  the  observed  and 
theoretical  data. 

In  our  case  the  observations  are  group  velocities  of 
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Love  and  Rayleigh  waves,  and  the  model  parameters  are 
shear  velocities  and  densities  for  a  horizontally  layered 
earth.  The  partial  derivatives  in  equation  1  are  generated 
numerically  by  the  method  described  by  Rodi  et  al.  (1975) . 

The  classic  least  squares  solution  to  equation  2  is 
by  minimizing  the  euclidean  length  of  AAX-AY  or 

||  AAX-AY  ||  =  (AAX-AY)T(AAX-AY)  =  £ 

Given  equation  2,  the  solution  is  (Hanson  and 
Lawson,  1975) : 

AAX  =  AY 

ATAAX  =  atay 

~  T  1  T 

AX  =  (A  A)'  A  AY 

where  AX  is  the  estimated  correction  vector. 

A  suitable  inverse  to  this  problem  can  be  obtained 
by  using  the  singular  value  decomposition  (SVD)  of  the 
matrix  A. 

A  =  UAVT 

where  the  columns  of  U  are  the  eigenvectors  associated 
with  the  colums  of  A,  the  rows  of  V  are  the  eigenvectors 
associated  with  the  rows  of  A,  and  A  is  a  diagonal  matrix 
of  non-zero  eigenvalues  of  the  matrix  A. 
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If  we  let 

K  =  (ATA)_1AT  (3a) 

then 

HAAX  =  HAY  (3b) 

and 

AX  =  HA  X  (3c) 

and 

AX  =  HAY 

Lsing  the  SVD  we  obtain  (Jackson,  1972) : 

AX  =  (AV~lUT)AY  (3d) 

H  =  AV" 1 UT 

Here,  H  is  known  as  the  pseudo-inverse  of  equation  2. 

The  matrix  product  HA=R  is  known  as  the  resolution 
matrix,  and  is  a  measure  of  the  uniqueness  of  the  solution 
(Jackson,  1972) .  The  rows  of  R  are  called  the  resolving 
kernals.  From  equation  3c,  it  can  be  seen  that  an  element 

A  ^  4  f  K 

of  AX,  AX^ ,  may  be  interpreted  as  convolving  the  x  row 
of  R  with  the  vector  AX.  Therefore,  AX^  can  be  considered 
the  weighted  sum  of  nearby  values  (Jackson,  1972) . 


I 

lo 
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Because  AaA  is  nearly  singular,  a  number  of  problems 
arise:  (1)  The  solution  vector  becomes  large  which  can 

cause  the  problem  to  leave  the  region  of  linearity;  (2) 
the  solution  oscillates  with  each  iteration.  By  looking 
at  the  pseudo- inverse  H  we  can  see  the  cause  of  this 
instability. 

H  =  VA“1UT 


where  the  matrix  A”1  equals 


1-1  = 


0  x2  0 


0  «...  y 

n  : 

.  T 

and  the  A^s  are  the  eigenvalues  of  the  matrix  A  A.  If 
T 

A  A  is  near  singular,  one  or  more  of  the  eigenvalues  will 
be  approaching  zero.  From  equation  3  it  can  be  seen  that 
a  small  eigenvalue  will  cause  a  large  change  in  one  or  more 


values  of  the  correction  vector  AX.  Similarily,  it  has 
been  shown  that  the  variance  of  the  model  parameters  are 
inversely  proportional  to  A  (Jackson,  1972)  .  Thus  small 
A’s  produce  large  standard  deviations. 

There  are  two  main  approaches  to  stabilizing  the 


inversion  process.  One  method  is  to  examine  the 
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eigenvalue  spectrum  of  the  matrix  A.  Small  eigenvalues 
are  removed  when  the  variance  becomes  too  large.  But 
removing  the  eigenvalues  degrades  the  resolution.  This 
trade-off  between  resolution,  and  variance  has  been 
discussed  by  a  number  of  authors  (Jackson,  1972;  Wiggins, 
1972;  and  Crosson,  1976) .  Braille  and  Keller  (1975)  used 
this  method  for  the  inversion  of  group  velocity  data. 

Inversion  stability  can  also  be  achieved  by  using 
Marquardt's  method  (Marquardt,  1963)  or  the  stochastic 
inverse  (Franklin,  1970) .  These  procedures  lead  to  the 
suppression  of  small  eigenvalues.  Franklin's  stochastic 
inverse  estimates  the  parameter  in  the  presence  of  noise. 
Equation  2  becomes: 

AAX  +  n  =  AY  (4) 

there  n  is  a  vector  of  observational  noise.  Franklin 
(1970)  shows  that  the  solution  to  equation  4  is  given  by: 

AX  =  WAT(AWAT  +  ED)“‘  AY  (5) 

where  W  is  the  covariance  matrix  of  the  parameter  and  ED 
is  the  covariance  matrix  of  the  observations. 

Marquardt's  method  considers  a  minimization  of  the 
functional 

(AAX  -  AY)T(EDTED) (AAX  -  AY)  -  AXT(WTW)  AX  (6) 
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with  the  solution  given  by: 

AX  =  (ATED_1  A  +  W  ~ 1 ) ~ 1 ATED" 1  AY  (7) 

kr 

.  T  i  T 

Letting  A  LD~  =  A  ,  equation  7  becomes: 

AX  =  (ATA  +  C2I)  ATAY  (3) 

where  c2  is  the  variance  of  the  parameter. 

By  looking  at  the  SVD  of  equation  8  we  can  see  how 
Marquardt's  method  tapers  the  eigenvalue  spectrum: 

AX  =  V| (A2  +  a 2 1 ) ~ 1 A | UTAY  (9) 

An  element  in  brackets  from  equation  9  equals:  (Crosson, 
1976) 


((A2  +  o2I)”'A|  .  =  X./U2.  +  a2) 

So  as  goes  to  0,  the  parameter  AX^  goes  to  zero.  A 
similar  analysis  for  the  stochastic  inverse  shows  that  it 
too  suppresses  small  eigenvalues. 

The  quantity  a2  can  also  be  thought  of  as  a  trade¬ 
off  parameter  between  resolution  and  variance  (Der  et 
al. ,  1970;  Wiggins,  1972;  Crosson,  1976).  In  Marquardt's 
method,  the  resolution  equals  (Crosson,  1976) : 

R  =  HA  =  (ATA  +  0 2 1 ) A^A  *  V|  (A  +  a2I)~'A2|VT  (10) 

For  the  stochastic  inverse 
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(11) 


R  =  I  (A2  +  O  2  I )  "  1  A  2  i  WT 

The  quantity  in  brackets  of  equations  10  and  11  is  a 
diagonal  matrix  whose  elements  have  the  form: 

Xi2/(X2  +  o2) 

If  o2  =  0,  then  R  =  I,  but  as  a2  goes  to  0  the  variance 
of  the  parameters  becomes  larger.  So  we  adjust  the  size 
of  a2  until  there  is  an  acceptable  trade-off  between 
standard  deviation  and  resolution. 

Two  other  elements  that  must  be  analyzed  in  the 
inversion  process  are  the  apriori  (SDX0)  and  postpriori 
(SOX)  standard  deviations  of  the  inversion  parameters: 

SDX0  =  DT(oWP“1)D 

SDX  = 

where  D  =  a  row  vector  of  a  delta  matrix  =  (0,0, . . . ,1, . . .0). 
This  analysis  will  give  a  measure  of  how  much  new  informa¬ 
tion  is  coming  from  the  data.  If  SDX0  and  SDX  are 
approximately  the  same,  we  know  that  very  little  informa¬ 
tion  is  coming  from  the  data. 

The  advantages  of  Marquardt’s  method  and  the  stochas¬ 
tic  inverse  is  that  the  SVD  of  A  is  not  explicitly  deter¬ 
mined  and  there  is  no  decision  to  be  made  on  the  rank  of 
T 

the  matrix  A  A. 

A  combination  of  these  two  methods  arc  used  in  our 
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inversion  process.  When  the  number  of  parameters  is 
greater  than  the  number  of  data,  the  stochastic  inverse 
is  used.  If  the  number  of  observations  is  greater  than 
the  number  of  parameters,  Marquardt's  method  is  used. 
This  procedure  leads  to  the  minimum  number  of 
computations . 


DISCUSSION 

Initial  models  for  the  inversion  were  based  on  the 
results  of  a  previous  surface  wave  study  done  by  Tung 
and  Teng  (1974).  Their  final  model  was  obtained  from 
Love  and  Rayleigh  wave  group  velocity  data.  Our 
starting  model  (Ml,  Table  5)  was  inverted  for  both  shear 
velocity  and  density. 

During  the  inversion,  large  instabilities  of  the 
density  parameter  were  found  in  the  crust.  This  was 
corrected  by  changing  the  starting  model  and  inverting 
for  compressional  wave  velocity.  In  the  derived  model 
(DM1)  a  low  velocity  layer  (LVL)  can  be  inferred  in  the 
sixth  layer,  but  the  validity  of  this  LVL  is  questionable 
because  of  the  insufficient  data  in  our  study.  Examin¬ 
ation  of  the  apriori  and  postpriori  standard  deviations 
for  the  shear  velocity  parameter  (Table  5)  shows  that 
there  is  only  a  small  amount  of  information  in  the  data 
applicable  to  depths  below  the  fifth  layer.  Thus  longer 
period  information  is  needed  to  confirm  velocities  beyond 


TABLE  5 


STARTING  MODEL 
Ml 


LAYER 


IICKNESS 

(KM) 

D-VELOCITY 

(KM/SEC) 

SHEAR  VELOCITY 
(KM/SEC) 

DENSITY 

10 

5.0 

3.0 

3.0 

10 

5.0 

3.1 

3.05 

20 

6.0 

3.5 

3.2 

20 

8.40 

4.6 

3.3 

25 

8,05 

4.55 

3.34 

50 

7.80 

4.2 

3.37 

50 

8.00 

4.3 

3.37 

FINAL  MODEL 
DM1 

SHEAR  VELOCITY 

(XM/SEC)  SDXo  SDX 


3.245 

.277 

.046 

3.193 

.277 

.076 

3.767 

.196 

.055 

4.498 

.196 

.106 

4.495 

.175 

.112 

4.283 

.124 

.105 

this  depth. 

A  new  initial  model  (M2)  using  thinner  layers  in 
the  crust  and  mantle  was  inverted  for  shear  wave  velocity, 
density  and  compress ion a 1  wave  velocity.  The  mantle 
shear  velocities  for  M2  were  also  modified,  as  were 
density  and  P  velocity  values  in  the  crust  (Table  6) . 

By  examining  the  relative  sizes  of  the  partial 
derivatives  of  different  parameters,  we  can  determine 
the  importance  of  each  parameter  in  the  inversion.  The 
partial  derivatives  of  group  velocity  with  respect  to 
P  wave  velocity,  shear  velocity  and  density  for  the 
fundamental  Rayleigh  wave  mode  are  plotted  in  Figure  11a, 
b,  and  c.  The  partial  derivatives  of  group  velocity  with 
respect  to  shear  velocity  and  density  are  plotted  in 
Figure  12a,  b  for  the  fundamental  Love  wave  mode.  Figures 
11  and  12  show  a  number  of  important  factors  that  affect 
this  study: 

(1)  The  partial  derivatives  for  P  wave  velocity  in 
the  upper  crust  are  not  insignificant  when  compared 
to  density  and  shear  velocity  and  therefore,  affect 
the  inversion. 

(2)  The  P  wave  partial  derivatives  in  the  lower 
crust  and  mantle  are  small,  and  are  not  significant 
in  the  inversion. 

(3)  The  partial  derivatives  with  respect  to  shear 
velocity  are  larger  than  those  for  density  and  P 
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TABLE  6 


STARTING  MODEL 
M2 


LAYER  P-VELOCITY  SHEAR  VELOCITY 


THICKNESS 

(KM/SEC) 

(KM/SEC) 

DENSITY 

5 

5 . 5 

3.0 

2.8 

5 

5.5 

3.0 

2.8 

5 

5.5 

3.1 

3.05 

5 

5.5 

3.1 

3.05 

10 

6.0 

3.5 

3.1 

10 

6.0 

3.5 

3.1 

10 

8.1 

4.6 

3.3 

10 

8.1 

4.6 

3.3 

10 

8.0 

4.4 

3.3 

15 

8.0 

4.4 

3.3 

25 

3.0 

4.4 

3.37 

25 

8.0 

4.4 

3.37 
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wave  velocity.  Therefore,  the  resolution  for  shear 
velocities  is  better  than  for  the  two  other 
parameters . 

The  starting  model  M2  and  the  derived  model  (SCSI) 
with  standard  deviation  error  bars  are  shown  in  Figure  13. 
The  fit  of  theoretical  group  velocities  for  SCSI  are 
shown  in  Figures  14  and  15.  The  averaging  kernals  are 
listed  in  Table  7,  where  only  the  elements  of  the  kernals 
corresponding  to  the  shear  parameters  are  included.  The 
other  elements  for  P  wave  velocity  and  density  are  small 
and  do  not  affect  the  structure.  The  width  of  these 
kernals  indicate  that  it  is  possible  to  resolve  layers 
to  a  thickness  of  5  km  in  the  upper  crust,  10  km  in  the 
middle  crust  and  20  km  in  the  lower  crust.  However, 
resolution  in  the  mantle  is  quite  poor.  All  the  standard 
deviations  were  found  to  be  less  than  .08  km/sec. 
Increasing  the  resolution  in  the  lower  crust  and  mantle 
(by  using  a  different  trade-off  parameter  o)  leads  to 
instability  in  the  shear  velocity  model,  and  unacceptably 
large  standard  deviations.  These  features  could  be 
improved  if  higher  mode  data  were  available. 

The  results  of  the  inversion  indicate  that  only  a 
few  layers  are  justified  to  describe  the  crust  with  the 
given  dispersion  information.  As  an  option  in  the  inver¬ 
sion  program,  averaging  kernals  can  be  made  to  resemble 


a  desired  box  car  vector.  An  example  would  be  if  we 
wished  to  average  the  first  two  parameters,  then  we  would 
want  the  kernal  to  resemble  a  box  car  function  such  as 
described  in  Table  8a.  In  the  inversion,  shear  wave 
parameters  from  layers  2,  3  and  4  were  averaged,  as  were 
parameters  from  layers  5  and  6,  and  layers  7  and  8, 

Table  8c.  The  actual  kernals  for  these  averaged  para¬ 
meters  are  shown  in  Table  8b.  The  averaged  values  and 
standard  deviations  are  plotted  in  Figure  16.  (Model 
ASCS1 ,  Table  9) .  A  comparison  with  Tung  and  Teng  (1974) 
results  indicates  here  that  crustal  shear  velocities 
are  higher  than  previously  reported  values.  And  the  lid 
velocity  in  the  mantle  is  found  to  be  slightly  less  than 
previously  determined.  A  LVL  at  60  km  is  consistent  with 
our  dispersion  data. 

A  list  of  the  apriori  and  postpriori  standard 
deviations  for  SCSI  is  given  in  Table  10.  There  is  almost 
no  new  information  that  can  be  derived  from  the  data  past 
70  km.  In  the  mantle  in  general,  the  amount  of  new 
information  from  the  data  is  quite  small.  But  the 
persistence  of  the  LVL  at  the  60-70  km  depth  might  indicate 
a  true  velocity  reversal  and  possibly  the  transition 
between  the  lithosphere  and  asthenosphere . 
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LAYER 

IICKNESS 

(KM) 

TABLE  9 

MODEL  ASCS  1 

SHEAR 

VELOCITY 

(KM/SEC) 

SDXq 

SDX 

5 

3.03 

.347 

.087 

15 

3.36 

.200 

.062 

20 

3.84 

.173 

.065 

20 

4.51 

.173 

.120 

10 

4.28 

.245 

.210 

15 

4.36 

.200 

.175 

25 

4.44 

.155 

.142 

25 

4.44 

.155 

.150 
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TABLE  10 


FINAL  MODEL 
SCP1 


SHEAR  VELOCITY 

SDXp 

SDX 

3.03 

.347 

.087 

3.43 

.347 

.169 

3.24 

.347 

.231 

3.4 

.347 

.234 

3.89 

.245 

.150 

3.78 

.245 

.169 

4.49 

.245 

.203 

4.53 

.245 

.210 

4.23 

.200 

.175 

4.36 

.155 

.142 

4. 

.155 

.150 

CONCLUSIONS 


By  comparing  the  size  of  the  P-wave  partial 
derivatives  with  density  and  shear  wave  partial  deriva¬ 
tives,  it  was  determined  that  compressional  wave  velocity 
cannot  be  neglected  when  inverting  surface  wave  data. 

The  resolution  for  shear  velocity  is  better  than  the 
resolution  for  density  and  p-wave  velocity. 

The  inversion  of  multi-layer  models  indicates  that 
a  3  layer,  40  km  crust  adequately  describes  the  South- 
China  subplate.  The  resolution  ranged  from  5  to  20  km  in 
the  crust  with  model  standard  deviations  of  less  than 
.08  km/sec.  The  shear  velocities  were  found  to  be  higher 
than  those  velocities  reported  by  Tung  (1974) . 

Upper  mantle  shear  velocities  were  found  to  be  lower 
than  previously  reported  values.  Resolution  in  the  mantle 
was  quite  poor  as  no  dispersion  information  is  obtained 
from  this  study  beyond  65  seconds.  The  apriori  and 
postpriori  standard  deviation  indicated  that  very  little 
new  information  was  being  obtained  from  the  data  about 
the  mantle.  But  because  of  the  persistence  of  a  LVL 
between  60  and  70  km,  it  is  presumed  that  transition 
between  the  lithosphere  and  asthenosphere  occurs  here. 
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